D s j(2317) meson production at RHIC 
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Production of D s j(2317) mesons in relativistic heavy ion collisions at RHIC is studied. Using 
the quark coalescence model, we first determine the initial number of D 3 j (2317) mesons produced 
during hadronization of created quark-gluon plasma. The predicted D s ,/(2317) abundance depends 
sensitively on the quark structure of the D s j(2317) meson. An order-of- magnitude larger yield is 
obtained for a conventional two-quark than for an exotic four-quark D s j (2317) meson. To include 
the hadronic effect on the Dsj (2317) meson yield, we have evaluated the absorption cross sections of 
the D s j (2317) meson by pion, rho, anti-kaon, and vector anti-kaon in a phenomenological hadronic 
model. Taking into consideration the absorption and production of D s ,/(2317) mesons during the 
hadronic stage of heavy ion collisions via a kinetic model, we find that the final yield of D s j (2317) 
mesons remains sensitive to its initial number produced from the quark-gluon plasma, providing 
thus the possibility of studying the quark structure of the D s j (2317) meson and its production 
mechanism in relativistic heavy ion collisions. 

PACS numbers: 25.75Nq, 12.39.Fe, 13.75.Lb,14,40.Aq 
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I. INTRODUCTION 



A narrow D s j(2317) meson was recently observed by 
the BABAR Collaboration [l[ in the inclusive Dfir in- 
variant mass distribution from e + e~ annihilation and 
confirmed by the Belle Collaboration in B meson decay 
0. This meson has the natural spin-parity J p = + 
and a mass below what was obtained from the QCD sum 
rule approach 0] and quark model calculations [ij for a 
normal two-quark state cs. The -D s j (2317) meson has 
thus been considered as a possible candidate for the ex- 
otic four-quark states that were studied in the bag model 
QCD sum rules Q and the non-relativistic poten- 
tial model §]. It is also possible that the D s j(2317) 
meson is simply a DK molecule or atom. Determination 
of the D s j(2317) meson width is limited by experimental 
resolutions to a value of less than 4.6 MeV/c 2 2]. The 
small width of -D s j(2317) meson is not surprising as its 
mass is below the threshold of DK system and can only 
decay into the kinematically allowed but isospin violated 
channel of D s tt state. Theoretically, the decay width of 
D s j — > D s n has been studied using the QCD sum rules, 
and its value varies from a few keV j| to a few tens 
keV [l(| depending on the assumed flavor state of four 
quarks or two quarks, respectively. A more phenomeno- 
logical approach based on the 3 Pq model [11| also gives 
a narrow width of about a few tens keV for a normal 
two-quark D s j (2317) meson. 

Studying the mechanism for D s j(2317) meson pro- 
duction in nuclear reactions is useful for understand- 
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ing its quark structure. In Ref.[12], a coupled-channel 
quark model has been used to study the production of 
a two-quark D s j (2317) meson and its radial excitations 
in hadronic reactions. Production of D s j(2317) in rel- 
ativistic heavy ion collisions has also been studied [l3l |. 
It was found that for a four-quark _D s j(2317) meson, a 
much large yield is obtained if one takes into account 
the diquark-diquark interactions in the produced quark- 
gluon plasma. Since the £) s j(2317) meson is not ex- 
pected to survive in the quark-gluon plasma, it is more 
likely to be produced at hadronization of the quark-gluon 
plasma either statistically or via quark coalescence. Its 
final abundance in a heavy ion collision depends, how- 
ever, also on its absorption and production probability 
in the subsequent hadronic matter. 

In the present paper, we study D s j (2317) meson pro- 
duction in central heavy ion collisions at the Relativistic 
Heavy Ion Collider (RHIC) in a kinetic model that starts 
from the final stage of the quark-gluon plasma, goes 
through a mixed phase of quark-gluon and hadronic mat- 
ters, and finally undergoes the hadronic expansion. The 
production of £> s j(2317) mesons from the quark-gluon 
plasma is modeled by the constituent quark coalescence 
model, which has been shown to describe reasonably not 
only the particle yields and their ratios [l4| but also 
their transverse momentum spectra and anisotropic flows 
ESESG3- The predicted number of £> s j(2317) mesons 
is found to depend on its quark structure, with the two- 
quark state giving an order-of-magnitude larger value 
than the four-quark state. The D s j(2317) meson can be 
absorbed and also produced in subsequent hadronic mat- 
ter via the reactions irD sJ «-> K*D(KD*), pD sJ «-> KD, 
KD sJ <-> pD(irD*), and K*D sJ ^ irD. The cross sec- 
tions for these reactions are evaluated in a phenomeno- 
logical hadronic model with coupling constants and form 
factors involving the D s j(2317) meson determined from 
the QCD sum rules. Taking into account the hadronic 
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effect in heavy ion collisions via a kinetic approach, we 
find that the final number of D s j(2317) mesons remains 
sensitive to its initial number produced from the quark- 
gluon plasma. Studying T> s j(2317) meson production in 
heavy ion collisions thus provides the possibility to study 
both its production mechanism and quark structure. 

This paper is organized as follows. In Section HH the 
dynamics of heavy ion collisions at RHIC is described. 
Production of Z? s j (2317) meson from the initial quark- 
gluon plasma via the quark coalescence model is dis- 
cussed in Section Mil The absorption cross sections of 
the D s j(2317) meson by hadrons such as the pion, rho, 
anti-kaon, and vector anti-kaon as well as their thermally 
averaged values are evaluated in Section Qv] Solving the 
rate equation based on a kinetic model, time evolution 
of the D s j (2317) meson abundance in heavy ion colli- 
sions is presented in Section [V] A summary is then given 
in Section IVII Details on the derivation of an approxi- 
mate analytical coalescence formula for D s j(2317) meson 
production from the quark-gluon plasma is given in Ap- 
pendix [3] while those on the QCD sum-rule approach to 
the determination of the form factor at the D s jDK ver- 
tex, which is needed in calculating the D s j (2317) meson 
absorption cross sections, is described in Appendix [B] 



smallpressure near phase transition and in hadronic mat- 
ter [2l| as well as to obtain a lifetime for the expand- 
ing matter comparable to that from the transport model 
[22| . Values of other parameters of the fire-cylinder are 
determined from fitting the measured transverse energy 
~ 788 GeV as well as the extracted freeze out tempera- 
ture Tp = 125 MeV and transverse flow velocity ~ 0.65c 
of midrapidity hadrons in central Au+Au collisions at 
^/snn = 200 GeV. Assuming that the fire-cylinder ex- 
pands isentropically, the fraction of hadronic matter dur- 
ing the mixed phase is found to increase approximately 
linearly, and the time dependence of the temperature of 
the fire-cylinder obtained in Ref.[l8j can be parameter- 
ized as 

T(r) =T C - (T H - T F ) ( 1 ^ L ) ^ , (2) 

\TF -Tff J 

where Th is the temperature of the hadronic matter at 
the end of the mixed phase and is thus the same as the 
critical temperature Tc for the quark-gluon plasma to 
hadronic matter transition. As in Ref. [l8j , we take Th = 
T c = 175 MeV. The freeze out temperature Tp = 125 
MeV then leads to a freeze out time Tp w 17.3 fm/c. 



II. HEAVY ION COLLISION DYNAMICS AT 
RHIC 



To model the dynamics of central relativistic heavy 
ion collisions after the end of the quark-gluon plasma 
phase, we use the schematic model of Ref. 18j based on 



the boost invariant picture of Bjorken 19( and an ac- 
celerated transverse expansion. Specifically, the volume 
of produced fire-cylinder in central Au+Au collisions at 
ytSNN — 200 GeV, which is the collision we are interested 
in, is taken to evolve with the proper time according to 

V{t) =n[R c +v c {T-T C ) + a/2{T-T C ) 2 } 2 tc, (1) 

where i?c = 8 fm and rc = 5 fm/c are final transverse 
and longitudinal sizes of the quark-gluon plasma, while 
vq = 0.4c is its transverse flow velocity at this time. 
The total transverse energy of quarks and gluons in the 
midrapidity (|y| < 0.5) is then about 1, 067 GeV if quarks 
and gluons are taken to be massive with m g = 500 MeV, 
rn u = rrid = m q = 300 MeV and m s = 475 MeV in or- 
der to take into account the nonperturbative effects of 
QCD near the critical temperature [2(j, and if the quark 
strangeness and baryon chemical potentials are taken to 
be /i s = and Lib = 10 MeV, respectively, to account 
for strangeness neutrality in the quark-gluon plasma and 
observed final antiproton to proton ratio of about 0.7 
at RHIC. The fire-cylinder then goes through a mixed 
phase of partonic and hadronic matters at a constant 
temperature Tc until th = 7.5 fm/c when its transverse 
radius and flow velocity are Rh ~ 9 fm and vh ~ 0.45c, 
respectively. This corresponds to a small transverse ac- 
celeration a = 0.02 c 2 /fm, which is chosen to reflect the 
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FIG. 1: (Color online) Time evolution of the abundance of 
pions, kaons, anti-kaons, and nucleons including contributions 
from decays of resonances (a) and the ratio of anti-baryon to 
baryon abundances (b) of mid-rapidity particles in central 
Au+Au collisions at ^/snn = 200 GeV. 

For normal hadrons such as pions, kaons, anti-kaons, 
and nucleons, they are taken to be in chemical equi- 
librium with baryon chemical potential lib = 30 MeV, 
charge chemical potential liq = MeV, and strangeness 
chemical potential n$ = 10 MeV. The nonzero strange 
chemical potential is needed to account for observed 
K~ / K + w 0.9 ratio in heavy ion collisions at RHIC. Ne- 
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glecting the time dependence of the chemical potentials, 
which have been shown to vary weakly with the temper- 
ature of an isentropically expanding matter in heavy ion 
collisions at RHIC [23[, time evolution of the abundance 
of pions, kaons, anti-kaons, and nucleons has been shown 
in Fig. 1 of Ref . [l8l| . Including also contributions from 
decays of resonances, time evolution of the abundance of 
these hadrons is shown in FigOJa). It is seen that the to- 
tal numbers of pions, kaons and anti-kaons do not change 
much during the hadronic stage while the nucleon num- 
ber decreases significantly as temperature drops. Their 
final numbers at freeze out are 926, 133, 113, and 47, 
respectively, and are comparable to those measured in 
experiments. In FigQJb), the time evolution of the ratio 
of anti-baryons to baryons is shown, and it changes from 
an initial value of 0.75 to a final value of 0.66, which is 
also comparable to the measured value of about 0.7. 



III. D sJ (2317) MESON PRODUCTION FROM 
THE QUARK-GLUON PLASMA 

A. The coalescence model 

In the coalescence model, the number of D s j (2317) 
mesons that are produced from the quark-gluon plasma 
is given by the product of a statistical factor gD B j which 
denotes the probability of combining cs or csqq quarks 
into a color neutral, spin 0, and isospin hadronic 
state and depends on whether the D s j (2317) meson is 
a two-quark or four-quark state, and the overlap of the 
quark phase-space distribution function f q (xi 1 pi) in the 
fire-cylinder with the Wigner distribution function 
of the D s j (2317) meson. The latter corresponds to 
the probability of converting the above partonic state 
into D s j (2317) and is dependent of the quark spatial 
wave functions in the D s j(2317) meson. Explicitly, the 
D s j(2317) number is expressed as 



N 



coal 



= 9D aJ 



nPi ■ da l d?Tp l 
fq{Xi,Pi) 



XfD Bj ( x l- x n;Pl-Pn), 



(3) 



with n — 2 or 4 for a two-quark or four-quark D s j (2317). 
Similar expressions have previously been used for study- 
ing the production of strange hadrons [24], charmed 
mesons [25j, and penta-quark baryons [18 1 from the 
quark-gluon plasma formed in relativistic heavy-ion col- 
lisions. 

We note that the coalescence model can be viewed as 
formation of bound states from interacting particles with 
energy mismatch balanced by other particles in the sys- 
tem. Neglecting such off-shell effects is reasonable if the 
binding energy is not large and / or if the production pro- 
cess is fast compared to the inverse of the energy mis- 
match. 

Since the normalized quark wave function of the 
-D s j(2317) meson in the color-spin-isospin space can be 



expressed as a linear combination of all possible or- 
thogonal flavor, color, and spin basis states, with co- 
efficients depending on the quark model used for the 
-D s j(2317) meson, gD aJ is simply given by the proba- 
bility of finding these quarks in any one of these color- 
spin-isospin basis states, i.e., go SJ — 1/3 2 x 1/2 2 = 1/36 
or 1/3 4 x 1/2 4 = 1/1296 for a two-quark or four-quark 
D s j(2317) meson, respectively. 

The da in Eq.© denotes an element of a space- like 
hypersurface ac at hadronization [26[ . In terms of the 
proper time r = (t 2 ~ z 2 ) 1 / 2 , the longitudinal momentum- 
energy rapidity y = ^ ln( E ^ ) z ) and space-time rapidity 

f]= \ lnCfz 2 )) the polar coordinates r and <f> in the trans- 
verse plane, the covariant volume element can be written 
as p ■ da — rrriT cosh(y — rj)rdrd(pdrj. 

For the phase-space distribution functions of quarks 
in the fire-cylinder, they are taken to be the same as in 
the description of the heavy ion collision dynamics, i.e., 
they are uniformly distributed in the transverse plane 
and their momentum distributions are relativistic Boltz- 
mannian in the transverse direction but uniform in ra- 
pidity along the longitudinal direction. Also, the Bjorkcn 
correlation of equal spatial 77 and momentum y rapidities 
is imposed, which is consistent with the small difference 
\y — 7/| < 0.5 seen in the transport model [22| . Explicitly, 
the quark momentum distribution per unit rapidity at 
T c is 



f q (v,r,y,PT) = g q S(v-y) 

-j(m T - p T • /3) - jJ,g 



x exp 



,(4) 



where g q — 6 is the color-spin degeneracy of a quark, 
(3 = (r/R)vc is the radial-dependent transverse flow ve- 
locity, 7 = (1 — /3 2 ) -1 / 2 , and the quark chemical potential 
Mg = Mh + Ms- The abundances of quarks at the end of the 
QGP phase in central Au+Au collisions at ^/snn = 200 
GeV described by the expanding fire-cylinder model of 
previous section are N u — Nd — 245, and N s ~ 149 at 
re , if we take into account the effect of gluons by convert- 
ing them into quarks according to the quark flavor com- 
position in the quark-gluon plasma as in Ref. [151]. For 
charm quarks, we assume that they are in thermal equi- 
librium in the quark-gluon plasma, which is supported 
by the large elliptic flow of electrons from charmed me- 
son decays that are observed in experiments [13, [28| and 
transport models [29l . [30j . Their number is, however, 
at present uncertain as it depends on whether charm 
quarks can be produced from the quark-gluon plasma. 
If we assume that the latter contribution is unimpor- 
tant, then the number of charm quarks N c produced in 
heavy ion collisions is simply given by the product of 
the charm quark number N^ N produced from an ini- 
tial hard nucleon-nucleon scattering and the total num- 
ber of binary collisions (about 960) in central Au+Au 
collisions. Using N* N from the PYTHIA program, we 
obtain iV c ~ 1.5. The value of N c increases to about 3 
and 7 if we use the N$ N from the PHENIX [U and 
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the STAR experiment [32j, respectively. In the following 
study, we will use N c = 3 for the calculation and discuss 
the sensitivity of our results to the change in the value for 
N c . For the charm quark mass, it is taken to be m c = 1.5 
GeV. 

For the Wigner distribution function of the Z? s j (2317) 
meson, instead of using a function of Lorentz invariant 
four-dimensional relative coordinates and momenta such 
as in Refs.pjl [Hj], we take it to be a function of three- 
dimensional relative coordinates and momenta for sim- 
plicity, based on the harmonic oscillator wave functions. 
Specifically, it is obtained from assuming that the wave 
functions of the quarks are those of a harmonic oscil- 
lator with an oscillator frequency u>. For a two-quark 
D s j(2317) meson with J n = + , its cs quarks are in the 
relative p-wave and its Wigner distribution function is 

thus m 



16 y 2 



16 

T 



x exp 



(5) 



In the above, we have used the usual definitions y 
x i — X2 and k = (m s pi — m c p 2 )/(rn c + m s ) for the 
relative coordinate and momentum of the two quarks, 
respectively. For the width parameter a, it can be related 
to the oscillator frequency uj by a = X/^lo) 1 / 2 with the 
reduced mass fi given by \i — m c m s / (m c + m s ). 

If the Z? s j (2317) meson is a four-quark meson, its four 
quarks are all in relative s- waves, which then leads to the 
Wigner distribution function 

P ) =8 3 exp (-E4-E • (6) 

\ i=l "i i=l / 

where the relative coordinates and momenta are 
related to the quark coordinates Xj and momenta by 
the Jacobian transformations: 



yi 



y2 



Y3 



xi - x 2 
2 / m r 



xi 



3 V w c + m s m c + m s 



-x 2 - x 3 



xi 



4 \ m c + m s + m q m c + m s + m 



-x 2 



m„ 



m c + m s + m q 



-x 3 - x 4 



(7) 



and 



ki = V2 
k 2 = 
k 3 = 



m s pi - m c p 2 
m c + m s 

3 m q (pi + p 2 ) - (m c + m s )p 3 

2 m c + m s + m q 

4 m q (pi + p 2 + p 3 ) - K + m s + m q )p 4 

3 m c + m s + 2m q 



It can be shown that the product of the Jacobians for 
the coordinate and momentum transformations is equal 
to one. 

The width parameter <Ji for the i-th relative coordinate 
in a four-quark D s j(2317) meson is again given by <7i — 
l/(/itiw) 1 / 2 with the reduced masses 



Mi 



/'2 



/'3 



2m c m s 
m c + m s ' 

3 m q (m c + m s ) 

2 m c + m s + m q ' 

4 m q {m c + m s + m g ) 

3 m c + m s + 2m q 



(9) 



We note that the reduced mass fi± of the cs quark pair 
in a four-quark D s j(2317) meson is a factor of two larger 
than that in a two-quark D s j (2317) meson due to dif- 
ferences in the definitions of the relative coordinate and 
momentum. 



B. Number of D s ./(2317) mesons produced from 
the quark-gluon plasma 

To evaluate the number of D s j(2317) mesons produced 
from the quark-gluon plasma requires information on the 
oscillator frequency to through the width parameter a 
in the D s j (2317) meson Wigner distribution function, 
which is related to the size of £> s j(2317). Since the latter 
is not known empirically, we choose the value of the oscil- 
lator frequency to fit instead the root-mean-square charge 
radius of the s-wave charmed Df(cs) meson. Taking its 
Wigner distribution function similar to Eq.® but with 
only one relative coordinate and momentum, we obtain 
the following mean-squared charge radius for the (cs) 
meson: 



( r D e )ch 



f<(xi-Y) 2 ) 

m 2 + 2m 2 
3(to c -I- m s ) 2 



2m 2 



2(to c + m s ) 2 



a 2 .(10) 



In the above, Y = (rn c xi + m s x 2 )/(m c + m s ) is the 
center-of-mass coordinate of cs quarks, and we have used 
the relation (y 2 ) = (3/2)<r 2 between the mean-square 
distance and the width parameter for two quarks in 
the relative s-wave as in the D+(cs) meson. Using the 
value (rfj ) c h w 0.124 fm 2 determined from the light- 
front quark model [34|, we find that a ~ 0.60 fm and 
huj re 300 MeV. 

For a two-quark D s j(2317) meson, whose quarks are 
in the relative p-wave, the relation between the mean- 
square distance of the two quarks and the width param- 
eter is (y 2 ) = (5/2)er 2 . Using above determined width 
parameter a, we obtain the following root-mean-square 
radius for a two-quark D s j(2317) meson: 



,1/2 

r / tWO 



1 (mg + m 2 ) 1 / 2 2 1/2 
y/2 m c + m 3 



5 



^5 (ml + m 2 s ) 1/2 



a « 0.53 fin. (11) 



For a four-quark D s j (2317) meson, its three size 
parameters are G\ = l/^/iiw) 1 / 2 w 0.42 fm, oa = 
l/^w) 1 / 2 w 0.58 fm, and ct 3 = l/^w) 1 / 2 « 0.6 fm. 
The resulting root-mean-square radius of a four-quark 
D s j(2317) meson is then 



r 2 \ 1/2 
D.j/four 



"3 (m 2 + m 2 )a 2 
4 (m c + TO S ) 2 
9 ((m c + rn s ) 2 + 2m 2 )(7| 



16 (m c 



m g ) 2 



1 ((m c + m s + m q ) 2 + 3m 2 )cr| 

2 (m c + to s + 2to 9 ) 2 
0.62 fm. 



1/2 



(12) 



which is somewhat larger than that of a two-quark 
D s j (2317) meson. 

The coalescence integral in Eq.([3]) can be evaluated 
analytically if we expand the hyperbolic functions to 
first order, neglect the transverse flow, and use non- 
relativistic momentum distributions for quarks. The first 
approximation is valid for \y\ < 0.5 considered in present 
study. Although the transverse flow affects strongly the 
transverse momentum spectrum of produced D s j(2317) 
mesons, it only has a small effect on its number. As 
shown in Appendix ^ these approximations lead to 
the following numbers of produced -D s j(2317) mesons 
from quark coalescence: ~ 1.9 x 10 ~ 2 for a two-quark 
-D s j(2317) meson and ~ 1.1 x 10~ 3 for a four-quark 
-D s j(2317) meson. These numbers are about a factor 
of two larger than those obtained from numerically eval- 
uating the coalescence integral using the Monte Carlo 
method of Ref . [l{| . which gives about 9.8 x 10 -3 and 
4.2 x 10 -4 per unit rapidity for the two-quark and four- 
quark D s j(2317) mesons, respectively, largely due to the 
use of the relativistic quark distribution functions. 

Since equilibrium thermal models have been success- 
fully employed in describing the experimental data for 
the yields and ratios of many hadrons in heavy ion col- 
lisions at RHIC [35l . [36| , it is of interest to compare the 
predicted number of D s j(23l7) mesons from the coales- 
cence model with that from the statistical model. In 
terms of the charm fugacity 7c and the strangeness chem- 
ical potential /is, this model gives the following number 
of produced £> s j(2317) mesons at hadronization: 



1C 



p^ddf, d 3 p 
, h (2tt) 3 E 
V Hlc e^slT H 



fo SJ {x,p) 



(27r) 2 j dm Tm T e 
Jo ( 25^1) « 5.2 x 10" 



(13) 



where fD aJ (x,p) is the thermal distribution function of 
D s j(2317) mesons, given by an expression similar to 



Eq.((4]) for quarks and Iq is the modified Bessel func- 
tion. In obtaining the numerical value in the last line 
of above equation, we have used Vh ~ 1,908 fm 3 , 
T H = 175 MeV, f3 H = 0.3c, ^ s = 10 MeV, and the charm 
fugacity 7c ~ 8.4. The latter ensures that the numbers 
of charmed hadrons produced statistically at hadroniza- 
tion is same as the number of charm quarks N c in the 
quark-gluon plasma. Specifically, we have No w 1.1, 
N D , fn 1.5, N Ds m 0.31, and N\ c w 0.11, giving a total 
of about 3 charmed hadrons. We note that the number 
of Z) s j(2317) mesons produced in the statistical model is 
independent of its quark structure, contrary to that in 
the coalescence model in which the yield for a two-quark 
D s j(2317) meson is about a factor of twenty larger than 
that for a four-quark one. 



IV. HADRONIC EFFECTS ON THE D sJ (2317) 
MESON 



A. D s j(2317) meson absorption cross sections by 
hadrons 
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FIG. 2: Born diagrams for D s ./(2317) absorption by ir, p, K, 
and K* mesons. 

The abundance of Z? s j(2317) mesons can change dur- 
ing the expansion of the hadronic matter as a result of 
absorption by pions, rho mesons, anti-kaons, and vector 
anti-kaons. Neglecting reactions with a D s meson in the 
final states, which are suppressed as a result of the pres- 
ence of the isospin violated vertex D s jD s tt, we have the 
following reactions: 



ttD s j — > KD*(K* D), pD sJ ^KD, 
KD 8 j — > pD(nD*), K*D sJ ^nD, 



(14) 



as shown in Fig[5] for the lowest-order Born diagrams. 
The cross sections for these reactions can be evaluated 
using the interaction Lagrangians 

CpKK = ig P KK{Kfdf,K - df.KfK) ■ p^, 
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£ p dd = -igpDDiDrd^D - d^DrD) ■ p^, 

Ck-Kt: = igK*KirK;7-(Kd^Tr-d^KTf)+K.c, 

Cd*d. = -ig D *i hr D*?-(Dff l i?-& , Di?)+-H.c., 

£d sJ dk = 9d sJ dkKDD s j. (15) 

In the above, r are Pauli matrices for isospin; and tt and 
p denote the pion and rho meson isospin triplet, respec- 
tively; and K = (K+,K°) T and K* = {K*+,K*°) T de- 
note the pseudoscalar and vector strange meson isospin 
doublet, respectively. The isospin doublet pseudoscalar 
D and vector D* mesons are defined in a similar way. 
For coupling constants, we use g pp = 2.52 from the 
vector dominance model (VDM) [33, HU, 9k*k-k = 3-25 
[39l | and gwDn — 6.3 from the decay widths of K* 
and D* , respectively, and g P KK = 3.25 from the SU(3) 
symmetry |4l| . For the coupling constant gKDD SJ , it has 
been studied in the QCD sum rules, and its value depends 
strongly on the quark structure of D s j(2317) meson. The 
predicted values are gD aJ DK = 9.2 GeV if the D s j(2317) 
meson is a two-quark state |42l| and gD s ,,DK = 3.15 GeV 
if it is a four-quark state [43| ■ 

The amplitudes for the reactions shown in Fig[2] are 
given by 



Mi 



M 2 = 



M 3a = 



1 



T ijgD 3 jDKgK*KTr 



g D s j DK g D* D-x 



t — m 2 K 



(2p 2 - Pi) 



M c 4 ! 



t — m? D 



(2p 2 - Pi)t,e 



4 i 



T i gD s .,DKg P KK 



-(2p 4 -p 2 ) M e 



M c 2 I 



K 



M zb = 
Mi 

M 6 



1 



T i gD s .,DKg P DD- 



-(P2 - 2^3)^6: 



M c 2 ) 



Ml(p 2 ^ 

M 3 (p 2 <-> 



-P3;P3 <-> 
-Pi)- 



Pa), 



(16) 



In the above, the matrix element rfj takes into account 
the isospin states of the particles in a reaction, with a 
denoting those of isospin triplet tt and p mesons, and i 
and j those of isospin doublet K, K* , D, and D* mesons. 
The momenta p\ and P2 are those of initial state particles 
while P3 and pi are those of final state particles on the 
left and right side of a diagram. The usual Mandelstam 
variables are given by s — (pi + p 2 ) 2 , t = (pi — P3) 2 , and 

u = (Pi - Pi) 2 - 

To obtain the full amplitudes, one needs in principle to 
carry out a coupled-channel calculation in order to avoid 
violation of unitarity. Such an approach is, however, be- 
yond the scope of present study. To prevent the artifi- 
cial growth of the tree-level amplitudes with the energy, 
we introduce instead form factors at interaction vertices, 
which are taken to have the form [4JI 



A 2 



A 2 



(17) 



where q , taken in the center of mass, is the squared 
three momentum transfer for t and u channels, or the 



squared three momentum of either the incoming or out- 
going particles for s channel. For the cutoff parameter 
A, we use A = 1.3 GeV for vertices involving an off-shell 
K meson and A = 3.7 GeV for those involving an off- 
shell D meson. These values are determined from the 
QCD sum-rule calculations given in Appendix [Bl for the 
D s jDK three-point functions. Although the calculations 
are only for a two-quark D s j(23l7) meson, we use them 
also for a four-quark Z) s ,/(2317) meson as well as for other 
vertices in the diagrams in Fig.[5J We expect this to be a 
reasonable assumption as a study of the V(3872) meson 
in the QCD sum rules has indicated that both the form 
and the cut-off of its form factor are not significantly dif- 
ferent between a two-quark and a four-quark V(3872) 
meson (45j. 

The isospin- and spin-averaged cross section is then 
given by 



1 



PS 



QA-ksNiNs Pi 



J dn\M n \ 2 F\ 



(18) 



where |A^ n | 2 denotes the squared amplitude obtained 
from summing over the isospins and spins of both ini- 
tial and final particles, with F denoting the appropri- 
ate form factors at interaction vertices. The factors 
Nj = (2/i + 1)(2J 2 + 1) and N s = (25 x + 1)(2S 2 + 1) in 
the denominator are due to averaging over the isospins I\ 
and I2 as well as the spins Si and S2 of initial particles. 
The three-momenta in the center of mass of initial and 
final particles are denoted by pi and pf, respectively. 
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FIG. 3: (color online) Cross sections for the absorption of 
a four-quark D 3 j (2317) meson by tt, p, K, and K* mesons 
via reaction ttDjj -* KD* (K* D), pD sJ -> KD, KD s j -> 
pD(nD*), and K*D sJ -> irD. 

In Fig. [31 we show the absorption cross sections of the 
D s j(2317) meson by tt, p, K, and K* as functions of 
the total center-of-mass energy s 1 / 2 above the threshold 
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energy s^ 2 of a reaction for the scenario that it is a four- 
quark state. Aside near the threshold of a reaction, where 
the cross section can be very large or small depending 
on whether the reaction is exothermic or endothermic, 
most cross sections are less than 1 mb except the reac- 
tion ttD s j — ► K*D, which has a peak value of about 2 
mb. If the D s j(2317) meson is a two-quark state, its ab- 
sorption cross sections are about a factor of nine larger 
than corresponding ones shown in Fig. [3] as the coupling 
constant gD 3j DK is about a factor of three larger for a 
two-quark D s j (2317) meson than for a four-quark one. 

The -D s j (2317) meson can also be produced in the 
hadronic matter by the inverse reactions KD*(DK*) — > 
w_D sJ , KD pD sJ , pD(wD*) -> KD sJ , and ttD 
K*D s j, with cross sections related to those of absorp- 
tion reactions via the detailed balance relations. 
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B. Thermally averaged D s ./(2317) meson 
absorption cross sections 



In the kinetic model to be used in the next section for 
studying D s j(22>\7) meson absorption and production in 
hadronic matter, the thermally averaged cross sections 
are needed. In terms of the thermal distribution func- 
tions /i(p) of D s j(2317) mesons and other hadrons, the 
thermally averaged cross section (Tabbed for the reaction 
ab — > cd is given by [46j 



J d 3 p a d 3 p b f a (p a )f b (p b )a ab 



cdVab 



J d 3 p a d 3 p b f a (p a )f b (p b ) 
= [AalKsia^alKsiab)}- 1 

dz[z 2 - (a a + a b ) 2 ][z 2 - (a a - a b ) 2 ] 
x Ki{z)a(s = z 2 T 2 ), (19) 

with oti = n%i/T, zq — max(a a + a bl a c + ad), K% being 
the modified Bessel function, and v ab denoting the rela- 
tive velocity of initial two interacting particles a and 6, 



FIG. 4: (color online) Thermally averaged cross sections for 
the absorption of a four-quark D s j(2317) meson by tt, p, K , 
and K* mesons via reaction nD s j — > KD*(K*D), pD s .j — > 
KD, KD SJ -> pD^D"), and K*D sJ -» tyD. 



V. TIME EVOLUTION OF THE D sJ (2317) 
MESON ABUNDANCE IN HADRONIC MATTER 



A. Rate equation for D s j (2317) meson production 
in heavy ion collisions 

In terms of thermally averaged cross sections and the 
densities of 7T, p, K, and K* mesons, the time evolution of 
the D s j(2317) meson abundance in the hadronic matter 
is determined by the kinetic equation 



dN DsJ ( T ) 
dr 



RQGp{r) + ^2(<TaD sJ ^bcVaD sJ )n a a) (t) 



njj) 



Nd sJ (t) 



(21) 



Vab 



V (Pa ■ Pb) 2 



E a E b 



(20) 



The thermally averaged absorption cross sections of 
the D s j(2317) meson by 7T, p, K, and K* as functions 
of the temperature of the hadronic matter are shown in 
Fig. 2]for the case that D s j(2317) is a four-quark meson. 
It is seen that the thermally averaged cross section of the 
dominant reaction ttD s j — > K*D has values less than 0.6 
mb. This value is again about a factor of nine larger if 
the .D s j(2317) is a two-quark meson. 



where n a (r) and (t) are, respectively, the equilibrium 
densities of light meson type a and charmed meson type c 
in the hadronic matter at proper time r when its temper- 
ature is T according to Eq. (J2]) , while Nq (t) is the equi- 
librium number of D s j (2317) mesons given by Eq. (|13|) 
using the temperature and flow velocity at proper time 
r. Since hadronization of the quark-gluon plasma takes 
a finite time of th — tc — 2.5 fm/c, Z? s j(2317) mesons 
are produced from the quark-gluon plasma in the mixed 
phase, with a rate proportional to the volume of the 
quark-gluon plasma. This is included in Eq. (|2ip through 
the term Rqgp (t) . Since the fraction of the quark-gluon 
plasma during the mixed phase decreases almost linearly 
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with the proper time, we can approximately write 



Rqgp{t) 



0, 

(0) 



t c < r < t h ; 
Otherwise. 



(22) 



In the above, N^ 1 is the total number of D s j(2317) 
mesons produced from the quark-gluon plasma. In the 
following calculations, it is obtained either from the co- 
alescence model by evaluating Eq.Q numerically using 
the Monte-Carlo method or from the statistical model 
using Eq. (fT3]) . In writing Eq. ([2~Tj) . we have assumed that 
the total number of charmed hadrons is conserved during 
the evolution of the hadronic matter as charms are not 
likely to be produced and destroyed in the hadronic mat- 
ter because of their small production and annihilation 
cross sections [43, [48|, |49j . 
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FIG. 5: Time evolution of the D s ./(2317) meson abundance 
in central Au+Au collisions at y / sjvjv=200 GeV for differ- 
ent initial numbers of D s j(2317) mesons produced from the 
quark-gluon plasma. 



B. D s j meson yield in relativistic heavy ion 
collisions 

In Fig[5l the abundance of Dsj (2317) mesons in cen- 
tral Au+Au collisions at ^/snn = 200 GeV is shown 
as a function of the proper time of the fire-cylinder. 
Since the initial number of D s j (2317) mesons produced 
from the quark-gluon plasma via quark coalescence is be- 
low the equilibrium number for both the two-quark and 
four-quark Z? s j(2317) mesons, its number increases dur- 
ing the hadronic evolution as shown by the dashed and 
solid lines, respectively. The final number of Z? s j (2317) 
mesons is about 3.0 x 10~ 2 if the £> s j(2317) meson is a 
two-quark state and is about 6.0 x 10 -3 if it is a four- 
quark meson. Although the ratio (~ 5) between the final 



numbers for the two and four-quark D s j (2317) mesons is 
smaller than that (~ 20) for initially produced Z? s j(2317) 
mesons, it is still appreciable. This result differs signif- 
icantly from the predictions of the statistical model. In 
this case, the D s j(2317) meson number decreases slightly 
to 3.8 x 10 -2 during the hadronic evolution if it is a 
two-quark meson as shown by the dash-dotted line, and 
it remains essentially unchanged during hadronic evolu- 
tion if it is a four-quark meson as shown by the dotted 
line. Since the final yield of D s j(2317) mesons in the 
coalescence model is much smaller for a four-quark state 
than for a two-quark state and also that from the sta- 
tistical model, studying Z? s j(2317) meson production in 
relativistic heavy ion collisions thus provides the possibil- 
ity of understanding not only its production mechanism 
but also its quark structure. 

Above results are obtained by assuming that there are 
three charm quarks in the quark-gluon plasma, based on 
the number of charm quarks measured by the PHENIX 
collaboration in p+p collisions. If this number is in- 
creased by a factor of two, which is closer to that expected 
from the STAR experiment on d+Au collisions, the final 
D s j(2317) meson numbers in both the coalescence and 
statistical models are increased by about a similar factor. 
A similar reduction factor is seen in the final D s j(2317) 
meson numbers in heavy ion collisions if the total charm 
quark number is reduced by a factor of two as given by 
the PYTHIA program for p+p collisions. 



VI. SUMMARY 

Using the quark coalescence model, we have predicted 
the yield of D s j (2317) mesons in central Au+Au colli- 
sions at RHIC. Contrary to the prediction of the statis- 
tical model, the initial number of D s j(2317) meson pro- 
duced at the end of the quark-gluon stage of heavy ion 
collisions depends sensitively on whether it is a two-quark 
or a four-quark meson, with the former giving an order 
of magnitude larger number than the latter. To take into 
account the effects of absorption and production during 
subsequent hadronic evolution, we have used a hadronic 
model to evaluate the cross sections for the absorption 
of the D s j(2317) meson by pion, rho meson, anti-kaon, 
and vector anti-kaon in the tree-level Born approxima- 
tion. With empirical masses and coupling constants as 
well as form factors from the QCD sum rules, we have 
found that all these cross sections are small, except the 
reaction irD s j — > K*D, which has a peak cross section 
of about 2 mb, if the D s ,/(2317) is a four-quark meson, 
but they are about ten times larger if it is a conventional 
two-quark meson. Including these reactions in a kinetic 
model based on a schematic hydro dynamic description of 
relativistic heavy ion collisions, we have studied the time 
evolution of the abundance of D s j(2317) mesons in these 
collisions. Our results show that the large difference in 
the initial numbers given by the quark coalescence model 
for the two-quark and four-quark D s j(2317) mesons re- 
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mains appreciable at freeze out. On the other hand, the 
-D s j(2317) number determined from the statistical model 
is essentially unchanged if the D s j (2317) is a four-quark 
meson and only change slightly if it is a two-quark me- 
son. Studying Z? s j(2317) production at RHIC and also at 
the forthcoming LHC thus offers the possibility to under- 
stand its quark structure and production mechanism. To 
achieve this goal requires, however, accurate information 
on the number of charm quarks produced during initial 
hard nucleon-nucleon collisions, as increasing or decreas- 
ing the initial quark number by a factor affects the final 
D s j(2317) meson number by a similar factor. Also, it 
is important to know if charm quarks can be produced 
from the quark-gluon plasma as this would affect the ini- 
tial D s j(2317) mesons produced from the quark-gluon 
plasma as well. 
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produced from quark coalescence: 



n Dsj = 9DsJ n n* n — j^y-^m — • 

(A3) 

In the above, we have made use of the fact the Wigner 
function of D s j(2317) meson is factorable in both the rel- 
ative coordinates and momenta of its constituent quarks. 
Because of— 0.5 < y — r\ < 0.5, the momentum space 
integral in Eq. (|A3|> reduces to a two-dimensional one, as 
the momentum integral in the z-direction gives simply 
one. 

If the D s j(23l7) meson is a p-wave two-quark meson, 
evaluating the integrals in Eq. (|A3[) with its Wigner func- 
tion given by Eq.© leads to the following number of 
-D s j(2317) mesons produced from quark coalescence: 



U s J 



1 



(47r) 3 / 2 M T C (T 5 



-TV ]V-- 
27 c s V c (l + 2^T c a 2 ) 2 

1.9 x 10~ 2 , 



(A4) 



For a four-quark D s j(2317) meson with its Wigner 
function given by Eq.©, the number of D s j (2317) 
mesons produced from quark coalescence is 



N 



f O LI 



1 



N c N s (N u Nu + N d N s ) 



1296 

f} 1 V c {l + 2,1^1 



1.1 x 10" 



(A5) 



APPENDIX A: APPROXIMATE EVALUATION 
OF THE COALESCENCE INTEGRAL 

For a hypersurface of constant proper time and a dis- 
tribution with Bjorken correlation between y and r), we 
can expand the hyperbolic function in the coalescence in- 
tegral to first order in y and rj if we consider D Si j(2317) 
meson production at midrapidity with \y\ < 0.5. In this 
case, the invariant phase-space factor in Eq.© can be 
approximated by 

Pi ■ dai— 1 ~ d 3 Xid 3 pi. (Al) 

Neglecting the transverse flow and treating quarks non- 
relativistically, we can then use the relation 

n r. 2 K 2 V 2 

V Pl = + V ' ( A2) 

^2m{T 2MT ^ 2uiT' y ' 

i—l i—1 

where M = X)"=i TOi ' ^° ex P ress the quark Boltzmann 
momentum distribution functions in terms of the total 
K and relative momenta. Using also the total and 
relative coordinates of the quarks, we obtain the fol- 
lowing expression for the number of D s _j(23l7) mesons 



APPENDIX B: THE D sJ DK FORM FACTOR 

In this appendix, we compute the D s jDK form factor 
using the QCD rum rules [5l| • In this approach, the 
short-range perturbative QCD is extended by the Wil- 
son's operator product expansion (OPE) of the correla- 
tors, which results in a series in powers of the squared 
momentum with Wilson coefficients. The convergence 
at low momentum is improved by using a Borel trans- 
form. The expansion involves universal quark and gluon 
condensates. Equating the quark-based calculation of a 
given correlator to the same correlator that is calculated 
using hadronic degrees of freedom via a dispersion rela- 
tion then provides the sum rules from which a hadronic 
quantity can be estimated. 

We shall use the three-point function to evaluate the 
D s jDK form factor by following the procedure suggested 
in Ref. [13] and further extended in [53[. This means 
that we shall calculate the correlators for an off-shell D 
meson and then for an off-shell K meson, requiring that 
the corresponding extrapolations to the respective poles 
lead to the same unique coupling constant. 

The three-point function associated with a D s jDK 
vertex with an off-shell D meson is given by 

r< D W) = Jd^xd^ioiTij^My^jom 
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(Bl) transform in both variables P 2 



where j' 5jU = sj^jsq, j D = iqj 5 c and j.d„j = cs are the 
interpolating fields for the K, D, and D s j, respectively 
with q, s and c being the light, strange and charm quark 
fields. Here we take the D s j to be a standard scalar 
quark-antiquark meson. 

The phcnomenological side of the vertex function, 
T^(p,p'), is obtained by the consideration of K and D 
states contribution to the matrix element in Eq. (|B1|) : 



m D s .,m 2 D F K f D f DsJ 



m c (p 2 



2 m 2 



K> 



9d sJ dk(t) , 
X (q 2 -m 2 D ) Pti 

+ higher resonances . (B2) 
In deriving Eq. (|B2|) . we have made use of 

(D sJ (p)\K(pi)D(q)) = g<g jDK (<l 2 ), (B3) 

where q — p' — p, and the decay constants Fk and fo 
and Jd sJ arc defined by the matrix elements 



(0\h»\K(p'))=ip'F K> 



>\3D\D(q)} = 



m 2 D f D 



in,;; 



and 



(Q\jD SJ \D sJ (p)) = m DsJ f DsJ 



(B4) 



(B5) 



(B6) 



The contribution of higher resonances and continuum in 
Eq. (|B2p will be taken into account as usual in the stan- 
dard form of Ref. j&ij , through the continuum thresholds 
So and uq for the D s j and K mesons, respectively. 

The QCD side, or theoretical side, of the vertex func- 
tion is evaluated by performing Wilson's operator prod- 
uct expansion of the operator in Eq. (|B1[) . Expressing 
r„ in terms of the invariant amplitudes, 



T ll {p,p l )^F l {p 2 1 p'\q 2 )p^ + F 2 {p 2 1 p'\q 2 )p^ (B7) 

we can write a double dispersion relation for each one of 
the invariant amplitudes, Fi, over the virtualities p 2 and 
p' 2 holding Q 2 = — q 2 fixed: 



= 



, p l (s,u,Q 2 ) 

ds 1 du l iTT 7i\ ' ( B8 ) 

o (s-p z ){u-p' ) 



where pi(s,u,Q 2 ) equals the double discontinuity of the 
amplitude Fi(p 2 ,p' ,Q 2 ) on the cuts m 2 < s < oo and 
< u < oo, which can be evaluated using Cutkosky's 
rules. Finally, in order to suppress the condensates of 
higher dimension and at the same time reduce the influ- 
ence of higher resonances, we perform a double Borel 



-P 



M 2 and 



P' A = —p' 2 — > M' 2 . Equating the two representations 
described above, we obtain the following sum rule in the 
structure p'^: 



= (Q 2 +m 2 D ) 



m r < ss > e 



-ml/M 2 



ds I duexp(~s/Al 2 )exp(~u/M' 2 ) 
o 



xf(s,t,u)9(u - u) 
where t = —Q 2 and 
f(s,t,u) 



(B9) 



T72 ( r 



2[A(a,u,*)]V 
+ [(2rn^ + 2m c m s 
+ s(t + u - s)][\(s,u,t)]~ 



2m r m, — s- 



t + u)(m 2 (s — t + u) 
(BIO) 



with X(s,u,t) = s 2 + u 2 + t 2 — 2su ~ 2st — 2tu, and 
u m ax — s + t -m 2 - st/m 2 . 

We use the same parameters as in Ref. [Hj]: m s — 
0.15 GeV, m c = 1.26 GeV, F K = 0.16 GeV, m D = 
1.865 GeV, m K = 0.498 GeV, m DsI = 2.317 GeV, f D = 
0.23 GeV, f DsJ = 0.225 GeV, (ss) = 0.8(qq), with 
(qq) = -(0.245) 3 GeV 3 . For the continuum thresholds 
we take so = (6.3±0.1) GeV 2 and uq = (rajf + Ait) 2 with 
Ait = 0.5 GeV. 

We also use the same Borel window as in Ref. [42j: 
10 GeV 2 < M 2 < 20 GeV 2 and work at a fixed ratio 
M' 2 /M 2 = 0.64/m|, s7 . We find a good Borel stability 
in this region of the Borel mass. Fixing M 2 — 15 GeV 2 , 
we show in Fig. 6 by the filled circles the momentum 
dependence of g D DK 

Since the present approach can not be used at small 
values of Q 2 , extracting the cjd sJ dk coupling from the 
form factor requires extrapolation of the curve to the 
mass of the off-shell meson D. In order to do this, we fit 
the QCD sum-rule results with an analytical expression. 
We have obtained a reasonable fit using a monopole form: 



a {D) (O 2 



92.4 



Q 2 + 14.1 ' 



(Bll) 



where the numbers are in units of GeV . This fit is also 
shown in Fig. 6 by the solid line. From Eq. (|Bll[ ) we get 

9d sJ dk = 9 { d^ dk (Q 2 = = 8 - 7 - 

To check the consistency of above fit, we also evaluate 
the form factor at the same vertex, but for an off-shell 
kaon. In this case, we have to evaluate the three-point 
function 



(B12) 
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a 

Q 
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Q 2 (GeV 2 ) 



+u(m c m s (s + 1 — u) + s(— s + 1 + it) 
+m 2 (-2u(s - t + u) + m c m s (s -t+ 3u))](B14) 



Using now uq = (ran + A u ) 2 with A u — 0.5GeV and 
M 2 = -r 5 M 2 , wc find that the results are also rather 

m D sJ 

stable as a function of the Borel mass. Fixing M 2 — 
15 GeV 2 , we show in Fig. 6 by the squares the QCD 

sum-rule results for 9^\dk^Q 2 )- A good fit of these 
results can be obtained using an exponential form: 



9 { nJ jDK (Q 2 ) = 7.98 e-^ 175 , 



(B15) 



FIG. 6: Momentum dependence of the D s jDK form factors. 
The solid and dashed lines give the parametrization of the 
QCDSR results for g ( £\ DK (Q 2 ) (circles) and g^ jDK (Q?) 
(squares), respectively. 



Proceeding in a similar way, we obtain the following sum 
rule: 



Q2+m « r ds r due -siM* e - u/ M> 



rn ~ J u 



4tt 2 
x g(s,t,u), 



where u m i n — rn 2 and 



(B13) 



g(s,t,u) 



[A(a,u,t)]3/ 1 



[m c (s — t + 3u) 



where 1.75 is in units of GeV 2 , as shown in Fig. 6 
by the dashed line. From Eq. (|B15[) we get Qd„jDK — 
^DsjDK^Q 2 ~ ~ m i") = 9.1, m excellent agreements 
with both the result obtained from g^ jriT< {Q 2 — —"m^) 
above and the result obtained in Ref. [42j for this cou- 
pling constant: Qd sJ dk — 9.2 ± 0.5. 

Considering the uncertainties in the continuum thresh- 
olds, and the difference in the values of the coupling 
extracted when the D meson or the kaon is off-shell, 
our result for the D s jDK coupling constant is thus 
9d sJ dk - 8.9 ± 0.9. 



^From the parameterizations in Eqs. (|Blip and (|B15[) . 
we can also get information about the cut-off (A) in the 
form factors. We see that the cut-off is much bigger when 
the D meson is off-shell (A « 3.7 GeV) than when the 
kaon is off-shell (A « 1.3 GeV), in agreement with the 
results obtained in refs. I521I53II. 
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